A distinct, high-affinity, alkaline phosphatase facilitates occupation of P-depleted environments by marine picocyanobacteria.

Marine picocyanobacteria of the genera Prochlorococcus and Synechococcus, the two most abundant phototrophs on Earth, thrive in oligotrophic oceanic regions. While it is well known that specific lineages are exquisitely adapted to prevailing in situ light and temperature regimes, much less is known of the molecular machinery required to facilitate occupancy of these low-nutrient environments. Here, we describe a hitherto unknown alkaline phosphatase, Psip1, that has a substantially higher affinity for phosphomonoesters than other well-known phosphatases like PhoA, PhoX, or PhoD and is restricted to clade III Synechococcus and a subset of high light I-adapted Prochlorococcus strains, suggesting niche specificity. We demonstrate that Psip1 has undergone convergent evolution with PhoX, requiring both iron and calcium for activity and likely possessing identical key residues around the active site, despite generally very low sequence homology. Interrogation of metagenomes and transcriptomes from TARA oceans and an Atlantic Meridional transect shows that psip1 is abundant and highly expressed in picocyanobacterial populations from the Mediterranean Sea and north Atlantic gyre, regions well recognized to be phosphorus (P)-deplete. Together, this identifies psip1 as an important oligotrophy-specific gene for P recycling in these organisms. Furthermore, psip1 is not restricted to picocyanobacteria and is abundant and highly transcribed in some α-proteobacteria and eukaryotic algae, suggesting that such a high-affinity phosphatase is important across the microbial taxonomic world to occupy low-P environments.

Complete genome sequence of Enterococcus faecalis phage EF_RCK.

A lytic bacteriophage EF_RCK infecting Enterococcus faecalis was isolated from a water sample collected in a raw cockle storage container at Taman Ria market, Sungai Petani, Malaysia. The phage has a 57,848-bp double-stranded DNA genome harboring 107 protein-encoding genes and shares 90.9% nucleotide similarity with Enterococcus phage EFKL (Saphexavirus genus).

Coordinated transcriptional response to environmental stress by a Synechococcus virus.

Viruses are a major control on populations of microbes. Often, their virulence is examined in controlled laboratory conditions. Yet, in nature, environmental conditions lead to changes in host physiology and fitness that may impart both costs and benefits on viral success. Phosphorus (P) is a major abiotic control on the marine cyanobacterium Synechococcus. Some viruses infecting Synechococcus have acquired, from their host, a gene encoding a P substrate binding protein (PstS), thought to improve virus replication under phosphate starvation. Yet, pstS is uncommon among cyanobacterial viruses. Thus, we asked how infections with viruses lacking PstS are affected by P scarcity. We show that the production of infectious virus particles of such viruses is reduced in low P conditions. However, this reduction in progeny is not caused by impaired phage genome replication, thought to be a major sink for cellular phosphate. Instead, transcriptomic analysis showed that under low P conditions, a PstS-lacking cyanophage increased the expression of a specific gene set that included mazG, hli2, and gp43 encoding a pyrophosphatase, a high-light inducible protein and DNA polymerase, respectively. Moreover, several of the upregulated genes were controlled by the host's phoBR two-component system. We hypothesize that recycling and polymerization of nucleotides liberates free phosphate and thus allows viral morphogenesis, albeit at lower rates than when phosphate is replete or when phages encode pstS. Altogether, our data show how phage genomes, lacking obvious P-stress-related genes, have evolved to exploit their host's environmental sensing mechanisms to coordinate their own gene expression in response to resource limitation.

Bacterial defense systems exhibit synergistic anti-phage activity.

Bacterial defense against phage predation involves diverse defense systems acting individually and concurrently, yet their interactions remain poorly understood. We investigated >100 defense systems in 42,925 bacterial genomes and identified numerous instances of their non-random co-occurrence and negative association. For several pairs of defense systems significantly co-occurring in Escherichia coli strains, we demonstrate synergistic anti-phage activity. Notably, Zorya II synergizes with Druantia III and ietAS defense systems, while tmn exhibits synergy with co-occurring systems Gabija, Septu I, and PrrC. For Gabija, tmn co-opts the sensory switch ATPase domain, enhancing anti-phage activity. Some defense system pairs that are negatively associated in E. coli show synergy and significantly co-occur in other taxa, demonstrating that bacterial immune repertoires are largely shaped by selection for resistance against host-specific phages rather than negative epistasis. Collectively, these findings demonstrate compatibility and synergy between defense systems, allowing bacteria to adopt flexible strategies for phage defense.

The long and short of it: benchmarking viromics using Illumina, Nanopore and PacBio sequencing technologies.

Viral metagenomics has fuelled a rapid change in our understanding of global viral diversity and ecology. Long-read sequencing and hybrid assembly approaches that combine long- and short-read technologies are now being widely implemented in bacterial genomics and metagenomics. However, the use of long-read sequencing to investigate viral communities is still in its infancy. While Nanopore and PacBio technologies have been applied to viral metagenomics, it is not known to what extent different technologies will impact the reconstruction of the viral community. Thus, we constructed a mock bacteriophage community of previously sequenced phage genomes and sequenced them using Illumina, Nanopore and PacBio sequencing technologies and tested a number of different assembly approaches. When using a single sequencing technology, Illumina assemblies were the best at recovering phage genomes. Nanopore- and PacBio-only assemblies performed poorly in comparison to Illumina in both genome recovery and error rates, which both varied with the assembler used. The best Nanopore assembly had errors that manifested as SNPs and INDELs at frequencies 41 and 157 % higher than found in Illumina only assemblies, respectively. While the best PacBio assemblies had SNPs at frequencies 12 and 78 % higher than found in Illumina-only assemblies, respectively. Despite high-read coverage, long-read-only assemblies recovered a maximum of one complete genome from any assembly, unless reads were down-sampled prior to assembly. Overall the best approach was assembly by a combination of Illumina and Nanopore reads, which reduced error rates to levels comparable with short-read-only assemblies. When using a single technology, Illumina only was the best approach. The differences in genome recovery and error rates between technology and assembler had downstream impacts on gene prediction, viral prediction, and subsequent estimates of diversity within a sample. These findings will provide a starting point for others in the choice of reads and assembly algorithms for the analysis of viromes.

Viral metagenomics reveals diverse virus-host interactions throughout the soil depth profile.

IMPORTANCE: Soil viruses can moderate the roles that their host microbes play in global carbon cycling. However, given that most studies investigate the surface layer (i.e., top 20 cm) of soil, the extent to which this occurs in subsurface soil (i.e., below 20 cm) is unknown. Here, we leveraged public sequencing data to investigate the interactions between viruses and their hosts at soil depth intervals, down to 115 cm. While most viruses were detected throughout the soil depth profile, their adaptation to host microbes varied. Nonetheless, we uncovered evidence for the potential of soil viruses to encourage their hosts to recycle plant-derived carbon in both surface and subsurface soils. This work reasons that our understanding of soil viral functions requires us to continue to dig deeper and compare viruses existing throughout soil ecosystems.

Biological Characterization and Genomic Analysis of Novel Phages DLDT_So2 and BHDT_So9 Against Pseudomonas solanacearum, an Infectious Agent in Tomato in Vietnam.

Tomato (Solanum lycopersicum L.) is an important grown vegetable in Vietnam. Bacterial wilt caused by Pseudomonas solanacearum has been considered to be an important disease resulting in a harvest loss up to 90% and significant economic loss to farmers. In this study, two bacteriophages DLDT_So2 and BHDT_So9 specific to P. solanacearum were isolated. Morphological analysis indicated that DLDT_So2 and BHDT_So9 had podovirus morphology and were classified into Autographiviridae family. The latent period and burst size of DLDT_So2 was found to be approximately 120 min and 20.0 ± 2.4 virions per infected cell. Meanwhile, the latent period of BHDT_So9 was 140 min with a burst size of 11.5 ± 2.8 virions per infected cell. Of the 23 bacterial strains tested, the phages infected 7/11 strains of P. solanacearum and none of the other bacteria tested were susceptible to the phages. Stability of the phages at different temperatures, pHs, solvents was also investigated. The genomes of DLDT_So2 and BHDT_So9 are 41,341 bp and 41,296 bp and long with a total GC content of 63%, contains 48 and 46 predicted protein-encoding CDSs. No virulence or antibiotic resistance genes were found in the genomes, suggesting they would be useful biocontrol agents against P. solanacearum. Classification of the phage using average nucleotide identity, phylogenetic analysis was also carried out. The two phages represented new species when they had overall average nucleotide identity of < 95%. This is first report of the isolation and characterization of P. solanacearum-specific phages from tomato farms in Vietnam. Supplementary Information: The online version contains supplementary material available at 10.1007/s12088-023-01090-9.

Novel Escherichia coli-Infecting Bacteriophages Isolated from Uganda That Target Human Clinical Isolates.

Background: The antimicrobial resistance catastrophe is a growing global health threat and predicted to be worse in developing countries. Phages for Global Health (PGH) is training scientists in these regions to isolate relevant therapeutic phages for pathogenic bacteria within their locality, and thus contributing to making phage technology universally available. Materials and Methods: During the inaugural PGH workshop in East Africa, samples from Ugandan municipal sewage facilities were collected and two novel Escherichia coli lytic phages were isolated and characterized. Results: The phages, UP19 (capsid diameter ∼100 nm, contractile tail ∼120/20 nm) and UP30 (capsid diameter ∼70 nm, noncontractile tail of ∼170/20 nm), lysed ∼82% and ∼36% of the 11 clinical isolates examined, respectively. The genomes of UP19 (171.402 kb, 282 CDS) and UP30 (49.834 kb, 75 CDS) closely match the genera Dhakavirus and Tunavirus, respectively. Conclusion: The phages isolated have therapeutic potential for further development against E. coli infections.

Novel ribotype/sequence type associations and diverse CRISPR-Cas systems in environmental Clostridioides difficile strains from northern Iraq.

The environment is a natural reservoir of Clostridioides difficile, and here, we aimed to isolate the pathogen from seven locations in northern Iraq. Four of the sites yielded thirty-one isolates (ten from soils, twenty-one from sediments), which together represent ribotypes (RTs) 001 (five), 010 (five), 011 (two), 035 (two), 091 (eight), and 604 (nine). Twenty-five of the isolates (∼81%) are non-toxigenic, while six (∼19%) encode the toxin A and B genes. The genomes of eleven selected isolates represent six sequence types (STs): ST-3 (two), ST-15 (one), ST-107 (five), ST-137 (one), ST-177 (one), and ST-181 (one). Five novel RT/ST associations: RT011/ST-137, RT035/ST-107, RT091/ST-107, RT604/ST-177, and RT604/ST-181 were identified, and the first three are linked to RTs previously uncharacterized by multilocus sequence typing (MLST). Nine of the genomes belong to Clade 1, and two are closely related to the cryptic C-I clade. Diverse multiple prophages and CRISPR-Cas systems (class 1 subtype I-B1 and class 2 type V CRISPR-Cas systems) with spacers identical to other C. difficile phages and plasmids were detected in the genomes. Our data show the broader diversity that exists within environmental C. difficile strains from a much less studied location and their potential role in the evolution and emergence of new strains.

Genomic analysis of Anderson typing phages of Salmonella Typhimrium: towards understanding the basis of bacteria-phage interaction.

The Anderson phage typing scheme has been successfully used worldwide for epidemiological surveillance of Salmonella enterica serovar Typhimurium. Although the scheme is being replaced by whole genome sequence subtyping methods, it can provide a valuable model system for study of phage-host interaction. The phage typing scheme distinguishes more than 300 definitive types of Salmonella Typhimurium based on their patterns of lysis to a unique collection of 30 specific Salmonella phages. In this study, we sequenced the genomes of 28 Anderson typing phages of Salmonella Typhimurium to begin to characterize the genetic determinants that are responsible for the differences in these phage type profiles. Genomic analysis of typing phages reveals that Anderson phages can be classified into three different groups, the P22-like, ES18-like and SETP3-like clusters. Most Anderson phages are short tailed P22-like viruses (genus Lederbergvirus); but phages STMP8 and STMP18 are very closely related to the lambdoid long tailed phage ES18, and phages STMP12 and STMP13 are related to the long noncontractile tailed, virulent phage SETP3. Most of these typing phages have complex genome relationships, but interestingly, two phage pairs STMP5 and STMP16 as well as STMP12 and STMP13 differ by a single nucleotide. The former affects a P22-like protein involved in DNA passage through the periplasm during its injection, and the latter affects a gene whose function is unknown. Using the Anderson phage typing scheme would provide insights into phage biology and the development of phage therapy for the treatment of antibiotic resistant bacterial infections.

The expendables: Bioarchaeological evidence for pauper apprentices in 19th century England and the health consequences of child labour.

Child labour is the most common form of child abuse in the world today, with almost half of child workers employed in hazardous industries. The large-scale employment of children during the rapid industrialisation of the late 18th and early 19th centuries in England is well documented. During this period, the removal of pauper children from workhouses in cities to work as apprentices in rural mills in the North of England was commonplace. Whilst the experiences of some of these children have been recorded historically, this study provides the first direct evidence of their lives through bioarchaeological analysis. The excavation of a rural churchyard cemetery in the village of Fewston, North Yorkshire, yielded the skeletal remains of 154 individuals, including an unusually large proportion of children aged between 8 to 20 years. A multi-method approach was undertaken, including osteological and palaeopathological examination, stable isotope and amelogenin peptide analysis. The bioarchaeological results were integrated with historical data regarding a local textile mill in operation during the 18th-19th centuries. The results for the children were compared to those obtained from contemporaneous individuals of known identity (from coffin plates) of comparable date. Most of the children exhibited distinctive 'non-local' isotope signatures and a diet low in animal protein when compared to the named local individuals. These children also showed severe growth delays and pathological lesions indicative of early life adversity, as well as respiratory disease, which is a known occupational hazard of mill work. This study has provided unique insights into the harrowing lives of these children; born into poverty and forced to work long hours in dangerous conditions. This analysis provides a stark testimony of the impacts of industrial labour on the health, growth and mortality risk of children, with implications for the present as well as our understanding of the past.

Genome Characterization of Enterococcus faecalis Bacteriophage EFKL.

We characterized the complete genome of the lytic Enterococcus faecalis phage EFKL, which was isolated from a sewage treatment plant in Kuala Lumpur, Malaysia. The phage, which was classified in the genus Saphexavirus, has a 58,343-bp double-stranded DNA genome containing 97 protein-encoding genes and shares 80.60% nucleotide similarity with Enterococcus phage EF653P5 and Enterococcus phage EF653P3.

Gauge your phage: benchmarking of bacteriophage identification tools in metagenomic sequencing data.

BACKGROUND: The prediction of bacteriophage sequences in metagenomic datasets has become a topic of considerable interest, leading to the development of many novel bioinformatic tools. A comparative analysis of ten state-of-the-art phage identification tools was performed to inform their usage in microbiome research. METHODS: Artificial contigs generated from complete RefSeq genomes representing phages, plasmids, and chromosomes, and a previously sequenced mock community containing four phage species, were used to evaluate the precision, recall, and F1 scores of the tools. We also generated a dataset of randomly shuffled sequences to quantify false-positive calls. In addition, a set of previously simulated viromes was used to assess diversity bias in each tool's output. RESULTS: VIBRANT and VirSorter2 achieved the highest F1 scores (0.93) in the RefSeq artificial contigs dataset, with several other tools also performing well. Kraken2 had the highest F1 score (0.86) in the mock community benchmark by a large margin (0.3 higher than DeepVirFinder in second place), mainly due to its high precision (0.96). Generally, k-mer-based tools performed better than reference similarity tools and gene-based methods. Several tools, most notably PPR-Meta, called a high number of false positives in the randomly shuffled sequences. When analysing the diversity of the genomes that each tool predicted from a virome set, most tools produced a viral genome set that had similar alpha- and beta-diversity patterns to the original population, with Seeker being a notable exception. CONCLUSIONS: This study provides key metrics used to assess performance of phage detection tools, offers a framework for further comparison of additional viral discovery tools, and discusses optimal strategies for using these tools. We highlight that the choice of tool for identification of phages in metagenomic datasets, as well as their parameters, can bias the results and provide pointers for different use case scenarios. We have also made our benchmarking dataset available for download in order to facilitate future comparisons of phage identification tools. Video Abstract.

Phenotypic Characterization and Comparative Genomic Analysis of Novel Salmonella Bacteriophages Isolated from a Tropical Rainforest.

Salmonella infections across the globe are becoming more challenging to control due to the emergence of multidrug-resistant (MDR) strains. Lytic phages may be suitable alternatives for treating these multidrug-resistant Salmonella infections. Most Salmonella phages to date were collected from human-impacted environments. To further explore the Salmonella phage space, and to potentially identify phages with novel characteristics, we characterized Salmonella-specific phages isolated from the Penang National Park, a conserved rainforest. Four phages with a broad lytic spectrum (kills >5 Salmonella serovars) were further characterized; they have isometric heads and cone-shaped tails, and genomes of ~39,900 bp, encoding 49 CDSs. As the genomes share a <95% sequence similarity to known genomes, the phages were classified as a new species within the genus Kayfunavirus. Interestingly, the phages displayed obvious differences in their lytic spectrum and pH stability, despite having a high sequence similarity (~99% ANI). Subsequent analysis revealed that the phages differed in the nucleotide sequence in the tail spike proteins, tail tubular proteins, and portal proteins, suggesting that the SNPs were responsible for their differing phenotypes. Our findings highlight the diversity of novel Salmonella bacteriophages from rainforest regions, which can be explored as an antimicrobial agent against MDR-Salmonella strains.

Staying below the Radar: Unraveling a New Family of Ubiquitous "Cryptic" Non-Tailed Temperate Vibriophages and Implications for Their Bacterial Hosts.

Bacteriophages are the most abundant biological entities in the oceans and play key roles in bacterial activity, diversity and evolution. While extensive research has been conducted on the role of tailed viruses (Class: Caudoviricetes), very little is known about the distribution and functions of the non-tailed viruses (Class: Tectiliviricetes). The recent discovery of the lytic Autolykiviridae family demonstrated the potential importance of this structural lineage, emphasizing the need for further exploration of the role of this group of marine viruses. Here, we report the novel family of temperate phages under the class of Tectiliviricetes, which we propose to name "Asemoviridae" with phage NO16 as a main representative. These phages are widely distributed across geographical regions and isolation sources and found inside the genomes of at least 30 species of Vibrio, in addition to the original V. anguillarum isolation host. Genomic analysis identified dif-like sites, suggesting that NO16 prophages recombine with the bacterial genome based on the XerCD site-specific recombination mechanism. The interactions between the NO16 phage and its V. anguillarum host were linked to cell density and phage-host ratio. High cell density and low phage predation levels were shown to favor the temperate over the lytic lifestyle for NO16 viruses, and their spontaneous induction rate was highly variable between different V. anguillarum lysogenic strains. NO16 prophages coexist with the V. anguillarum host in a mutualistic interaction by rendering fitness properties to the host, such as increased virulence and biofilm formation through lysogenic conversion, likely contributing to their global distribution.

Abolishment of morphology-based taxa and change to binomial species names: 2022 taxonomy update of the ICTV bacterial viruses subcommittee.

This article summarises the activities of the Bacterial Viruses Subcommittee of the International Committee on Taxonomy of Viruses for the period of March 2021-March 2022. We provide an overview of the new taxa proposed in 2021, approved by the Executive Committee, and ratified by vote in 2022. Significant changes to the taxonomy of bacterial viruses were introduced: the paraphyletic morphological families Podoviridae, Siphoviridae, and Myoviridae as well as the order Caudovirales were abolished, and a binomial system of nomenclature for species was established. In addition, one order, 22 families, 30 subfamilies, 321 genera, and 862 species were newly created, promoted, or moved.

Functional ecology of bacteriophages in the environment.

Bacteriophages are as ubiquitous as their bacterial hosts and often more abundant. Understanding how bacteriophages control their bacterial host populations requires a number of different approaches. Bacteriophages can control bacterial populations through lysis, drive evolution of bacterial immunity systems through infection, provide a conduit for horizontal gene transfer and alter host metabolism by carriage of auxiliary metabolic genes. Understanding and quantifying how bacteriophages drive these processes, requires both technological developments to take measurements in situ, and laboratory-based studies to understand mechanisms. Technological advances have allowed quantification of the number of infected cells in situ, revealing far-lower levels than expected. Understanding how observations in laboratory conditions relate to what occurs in the environment, and experimental confirmation of the predicted function of phage genes from observations in environmental omics data, remains challenging.

Crop management shapes the diversity and activity of DNA and RNA viruses in the rhizosphere.

BACKGROUND: The rhizosphere is a hotspot for microbial activity and contributes to ecosystem services including plant health and biogeochemical cycling. The activity of microbial viruses, and their influence on plant-microbe interactions in the rhizosphere, remains undetermined. Given the impact of viruses on the ecology and evolution of their host communities, determining how soil viruses influence microbiome dynamics is crucial to build a holistic understanding of rhizosphere functions. RESULTS: Here, we aimed to investigate the influence of crop management on the composition and activity of bulk soil, rhizosphere soil, and root viral communities. We combined viromics, metagenomics, and metatranscriptomics on soil samples collected from a 3-year crop rotation field trial of oilseed rape (Brassica napus L.). By recovering 1059 dsDNA viral populations and 16,541 ssRNA bacteriophage populations, we expanded the number of underexplored Leviviricetes genomes by > 5 times. Through detection of viral activity in metatranscriptomes, we uncovered evidence of "Kill-the-Winner" dynamics, implicating soil bacteriophages in driving bacterial community succession. Moreover, we found the activity of viruses increased with proximity to crop roots, and identified that soil viruses may influence plant-microbe interactions through the reprogramming of bacterial host metabolism. We have provided the first evidence of crop rotation-driven impacts on soil microbial communities extending to viruses. To this aim, we present the novel principal of "viral priming," which describes how the consecutive growth of the same crop species primes viral activity in the rhizosphere through local adaptation. CONCLUSIONS: Overall, we reveal unprecedented spatial and temporal diversity in viral community composition and activity across root, rhizosphere soil, and bulk soil compartments. Our work demonstrates that the roles of soil viruses need greater consideration to exploit the rhizosphere microbiome for food security, food safety, and environmental sustainability. Video Abstract.

Antimicrobial resistance in dairy slurry tanks: A critical point for measurement and control.

Waste from dairy production is one of the largest sources of contamination from antimicrobial resistant bacteria (ARB) and genes (ARGs) in many parts of the world. However, studies to date do not provide necessary evidence to inform antimicrobial resistance (AMR) countermeasures. We undertook a detailed, interdisciplinary, longitudinal analysis of dairy slurry waste. The slurry contained a population of ARB and ARGs, with resistances to current, historical and never-used on-farm antibiotics; resistances were associated with Gram-negative and Gram-positive bacteria and mobile elements (ISEcp1, Tn916, Tn21-family transposons). Modelling and experimental work suggested that these populations are in dynamic equilibrium, with microbial death balanced by fresh input. Consequently, storing slurry without further waste input for at least 60 days was predicted to reduce ARB spread onto land, with > 99 % reduction in cephalosporin resistant Escherichia coli. The model also indicated that for farms with low antibiotic use, further reductions are unlikely to reduce AMR further. We conclude that the slurry tank is a critical point for measurement and control of AMR, and that actions to limit the spread of AMR from dairy waste should combine responsible antibiotic use, including low total quantity, avoidance of human critical antibiotics, and choosing antibiotics with shorter half-lives, coupled with appropriate slurry storage.

Elucidating the picocyanobacteria salinity divide through ecogenomics of new freshwater isolates.

BACKGROUND: Cyanobacteria are the major prokaryotic primary producers occupying a range of aquatic habitats worldwide that differ in levels of salinity, making them a group of interest to study one of the major unresolved conundrums in aquatic microbiology which is what distinguishes a marine microbe from a freshwater one? We address this question using ecogenomics of a group of picocyanobacteria (cluster 5) that have recently evolved to inhabit geographically disparate salinity niches. Our analysis is made possible by the sequencing of 58 new genomes from freshwater representatives of this group that are presented here, representing a 6-fold increase in the available genomic data. RESULTS: Overall, freshwater strains had larger genomes (≈2.9 Mb) and %GC content (≈64%) compared to brackish (2.69 Mb and 64%) and marine (2.5 Mb and 58.5%) isolates. Genomic novelties/differences across the salinity divide highlighted acidic proteomes and specific salt adaptation pathways in marine isolates (e.g., osmolytes/compatible solutes - glycine betaine/ggp/gpg/gmg clusters and glycerolipids glpK/glpA), while freshwater strains possessed distinct ion/potassium channels, permeases (aquaporin Z), fatty acid desaturases, and more neutral/basic proteomes. Sulfur, nitrogen, phosphorus, carbon (photosynthesis), or stress tolerance metabolism while showing distinct genomic footprints between habitats, e.g., different types of transporters, did not obviously translate into major functionality differences between environments. Brackish microbes show a mixture of marine (salt adaptation pathways) and freshwater features, highlighting their transitional nature. CONCLUSIONS: The plethora of freshwater isolates provided here, in terms of trophic status preference and genetic diversity, exemplifies their ability to colonize ecologically diverse waters across the globe. Moreover, a trend towards larger and more flexible/adaptive genomes in freshwater picocyanobacteria may hint at a wider number of ecological niches in this environment compared to the relatively homogeneous marine system.